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ABSTRACT 
A  method  for  nonlinear  system  design,  based  on  the  generalized 
Popov  stability  criterion,  is  developed.    It  is  shown  that,  although  the 
design  method  is  applicable  for  all  values  of  q  in  the  Popov  theory,  it  is 
particularly  useful  when  q  is  non-zero,  since  it  greatly  simplifies  the 
design  effort  for  this  case.    All  designs  are  accomplished  in  conjunction 
with  the  Nyquist  and  modified  Nyquist  loci.    Basic  to  the  design  proce- 
dure is  the  development  and  utilization  of  the  modified  frequency-response 
polar  loci  for  lag  and  lead  compensation  networks.    Three  examples,  one 
with  digital  simulation,  are  included  to  illustrate  the  procedure  described. 
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I  .     INTRODUCTION 

The  now  well-known  Popov  criterion  has  provided  a  new  approach 
to  the  idea  of  absolute  stability  in  the  dynamics  of  automatic  control 
systems.    As  a  result,  the  large  amount  of  subsequent  research  has 
contributed  much  to  the  development  of  modern  stability  theory. 

The  importance  of  stability  theory  is  quite  evident  when  one  realizes 
that  all  the  current  conventional  design  techniques  for  linear,  time-invar- 
iant systems  are  directly  or  indirectly  derived  from  stability  theory.    This 
is  reflected  in  the  fact  that  the  basic  engineering  design  approach  is  to 
determine  how  an  unstable  system  can  be  altered  to  make  it  stable. 
When  the  system  is  nonlinear    or  time-varying,  the  design  problem 
becomes  much  more  complex.    However,  considerable  progress  has  been 
made  in  this  area  in  recent  years  . 

This  paper  considers  systems  of  the  standard  form  of  Figure  1.    The 
linear  element  may  contain  time  delays  or  distributed  parameters.    The 
nonlinear  element  may  be  time-varying  and  may  contain  hysteresis.    The 
only  restriction  on  the  nonlinear  element  is  that  its  input-output  character- 
istic must  lie  within  a  sector  bounded  by  two  straight  lines  passing 
through  the  origin.     There  is  no  loss  of  generality  if  one  line  corresponds 
to  the  x-axis  since  this  can  be  accomplished  through  a  simple  transfor- 
mation.    Hence: 

0   <  -^   £  k         for  all  t.  (1) 
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This  class  of  systems  may  be  treated,  in  some  cases,  by  phase- 
plane  methods  or  by  approximate  methods  such  as  the  single  or  dual- 
input  describing  function.    Some  exact  methods  are  also  historically 
available. 

The  solution  of  the  absolute  stability  problem  was  first  formulated 
in  1944  by  Lur'e  and  Postnikov  [1]  for  the  case  when    u  =  f(e)      is  located 
completely  in  the  first  and  third  quadrants     (k  =  °°).     Also,  in  their  paper, 
the  utilization  of  a  Lyapunov  function  of  the  type  "a  quadratic  form  plus 
an  integral  involving  the  nonlinearity"  was  first  proposed.    In  further 
works,  Lur'e  developed  a  method  which  permits  one  to  write  down  a 
finite  system  of  quadratic  equations  directly  from  the  given  differential 
systems  .    He  showed  that  if  this  system  of  equations  ,  which  he  called 
resolving  equations,  has  real  roots,  then  this  is  a  sufficient  condition 
for  the  absolute  stability  of  the  system  of  Figure  1. 

Yakubovitch  [l  ]  further  developed  this  approach  by  showing  that  the 
system  of  resolving  equations  can  be  reduced  to  a  single  algebraic 
equation.     He  then  proved  that  the  existence  of  two  real  roots  for  this 
equation  is  a  sufficient  condition  for  absolute  stability. 

Most  of  the  subsequent  research  was  centered  around  the  application 
of  Lyapunov1  s  method  for  determining  system  stability.     However,  up 
until  this  time,  all  the  exact  methods  suffered  from  the  disadvantage 
that  even  the  solution  of  the  simple  problems  became  so  formidable 
computationally  that  the  methods  were  generally  unusable  as  engineering 
tools . 


Then,  in  1959,  Popov  proposed  an  entirely  different  approach  to  the 
problem  of  absolute  stability.     He  expressed  his  sufficient  conditions 
in  terms  of  the  frequency  response  of  the  linear  portion  of  the  system, 
thus  giving  his  criterion  a  very  simple  form  that  is  convenient  for  graph- 
ical application.    In  fact,  the  Popov  criterion  is  so  powerful  that  all  the 
results  described  above  and  connected  with  the  Lyapunov  function,  con- 
sisting of  a  quadratic  form  plus  an  integral  of  the  nonlinearity ,  are 
included  in  the  criterion. 

In  the  Popov  method  only  the  Nyquist  diagram  of  the  linear  element 
and  the  inequality  constraints  of  the  nonlinear  element  are  needed. 
Hence  the  criterion  is  readily  applied  to  higher-order  systems  and  those 
containing  time  delays  or  distributed  parameters  .    Since  it  is  basic  to 
the  design  method  developed  in  this  paper  a  brief  development  of  the 
Popov  theorem  follows.    For  a  proof  of  the  theorem,  see  [l]  . 

In  the  most  simple  case,  the  problem  of  absolute  stability  of  dynamic 
processes  consists  of  finding  when  the  equilibrium  state  of  the  general 
system 

x  =  Ax  +  bu 

u  =  f(e)  (2) 


where  x  is  a  vector,  A  is  a  square  matrix,  b  is  a  column  vector,  c'  is  a 
row  vector,  and  u  is  a  scalar,  is  on  the  whole  asymptotically  stable  in 
the  Lyapunov  sense.    The  formulation  (2)  is  a  system  of  differential 
equations  differing  from  a  linear  homogeneous  set  with  constant 
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coefficients  only  by  a  nonlinear  function  whose  argument  is  any  one  or  a 
linear  combination  of  the  state  variables. 

A,  b,  and  c  are  real,  time-invariant  matrices  in  which  some  of  the 
components  may  be  zero.    The  nonlinear  characteristic,     f(e)   ,    is  an 
arbitrary,  single-valued,  piecewise  continuous  real  function,  defined 
for  all  real  values  of    e    and  satisfying  the  conditions 

f(0)  =  0  (3) 

e 
which  are  interpreted  geometrically  in  Figure  2  . 

The  problem  of  establishing  the  absolute  stability  of  the  system  (2) 

requires  the  following  definitions  .     A  system  (2)  such  that  all  the  roots 

of  its  linear  portion  characteristic  equation  lie  in  the  left  half-plane  will 

be  the  so-called  principal  case  of  (2).     However,  a  case  when  some  of 

the  roots  of  the  characteristic  equation  are  on  the  imaginary  axis  and  the 

rest  are  in  the  left  half-plane  will  be  called  a  particular  case  of  (2).    In 

(2),  assume  that  f(e)  =  he  with  0  ^  h    <  k  ,    for  the  principal  case  and 

with    0    <  h   <  k    for  the  particular  cases.    Then  (2)  becomes  the  linear 

system 

x  =  Ax  +  hbe 

(4) 

c  =  c's   „ 

Now,  Popov  expressed  his  results  in  the  frequency  domain.     In  addition 
to  the  usual  open-loop  frequency  response,  G(jw)   (see  Figure  1),  he  used 
a  modified  frequency  response,  G*(jto),  which  is  defined  as  follows 
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u 


u  =  f  (e) 


Fig.  2.    The  input-output  characteristic  of  the  nonlinear 
element. 
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Re  G*(jw)  =  Re  G(joo)  =X 

(5) 
Im  G*(jw)  =  03  Im  G(jw)  =  Y  . 

From  linear  theory,  for  the  system  (4),  obtained  from  the  principal 
case  of  (2),  to  be  stable  for  any    h    in  the  interval    0  ^  h    <  k,    it  is 
necessary  and  sufficient  that  the  locus  of  G(jw) ,    and  hence  also  the 
locus  of  G*(joO,     never  intersect  the  real  axis  from  negative  infinity  up 
to  and  including  the  point  (-1/k)  .    Also  for  the  system  (4),  corresponding 
to  particular  cases  of  (2),  to  be  stable  with  0  <  h    £   k     ,    it  is  necessary 
and  sufficient  that  the  locus  of    G(joo)   (or  of  G*(jo>)  )  never  intersect  the 
forbidden  zone  (the  interval  of  the  negative  real  axis  from  negative 
infinity  to  -  l/k)  and  that  the  system  (4)  be  stable-in-the-limit,  that  is, 
that  the  system  (4)  be  stable  for  sufficiently  small  h    >  0     (that  is, 
h   =   c  >  0    where  c   is  arbitrarily  small) . 

Consider  now  the  original  problem  of  establishing  the  absolute 
stability  of  the  system  (2),  that  is,  when    f(e)    is  an  arbitrary  (not 
necessarily  linear)  characteristic  contained  in  the  sector  (0,k)  for  the 
principal  case,  or  in  the  sector    (c   ,    k)    for  the  particular  cases. 
Additional  conditions  must  be  imposed  on  the  frequency  response  in 
order  to  guarantee  stability  for  the  system  (2),  not  only  for  linear  charac- 
teristics   f(e)  =  he  ,    but  for  any  nonlinear  characteristics    f (e)    contained 
in  the  given  sector.     Popov  stated  these  additional  requirements  in  the 
following  famous  theorem  [l  ]  . 

Theorem.    For  the  system  (2)  to  be  absolutely  stable  in  the  sector 
(0,k)  for  the  principal  case,  and  in  the  sector  (c  ,  k)  for  the  particular 
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cases  (where  e>  0  is  an  arbitrarily  small  number),  it  is  sufficient  that 
there  exist  a  finite  real  number  q  such  that  for  all  a.'  ^  0  the  following 
inequality  is  satisfied 

Re(l  +  jqw)G(jO>)  +  l/k>  0  (6) 

and,  additionally  for  the  particular  cases,  that  the  conditions  for 
stability-in-the-limit  be  satisfied. 

In  order  to  give  the  criterion  a  geometric  interpretation  let 

G*(ja>)  =X  + jY  (7) 

Then 

Re(l  +  jqa>)G(jw)  =  ReG(jw)  -  qwImG(ju>)  =  X-qY. 
Hence  (6)  can  be  written  as 

X-qY  +  1/k  >  0    for  all  u>  *  0  .  (9) 

But  the  equation 

X-qY+l/k=0  (10) 

is  the  equation  of  a  straight  line  (the  Popov  line)  with  slope  1/q  which 
passes  through  the  point  -1/k  on  the  real  axis.    Hence  it  is  possible  to 
give  a  geometric  interpretation  to  the  Popov  theorem. 

For  the  system  (2)  to  be  absolutely  stable  in  the  sector  (0,k)  for  the 
principal  case,  or  in  the  sector  (c  ,  k)   for  the  particular  cases,  it  is 
sufficient  that  there  exists  in  the  G*(ju>)  plane  a  straight  line,  passing 
through  the  point  -1/k  on  the  real  axis,  such  that  the  modified  frequency 
response    G*(jw)  lies  strictly  to  the  right  of  it.    Also,  for  the  particular 
cases,  the  conditions  for  stability-in-the-limit  must  be  satisfied.    Some 
examples  are  shown  in  Figure  (3). 
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Attempts  to  extend  the  Popov  theorem  have  met  with  considerable 
success.    Hsu  and  Meyer  [2]  have  presented  a  generalized  Popov  theorem 
which  includes  the  recent  extensions  . 

Theorem.    Consider  the  basic  system  of  Figure  1  with  the  linear 
element  being  output  stable.    For  the  system  to  be  both  absolutely 
control  and  output  asymptotic  for  (u/e)  c  (0/k) ,  it  is  sufficient  that  there 
exists  a  real  number  q  such  that  for  all  real  a)  ^  0   and  an  arbitrarily 
small    6  >    0,  the  following  inequality  is  satisfied: 

Re  (  (1  +  ja)q)G(ja))  )  +  1/k  *  6  >  0  (11) 

The  restrictions  on  qand  k  depending  on  the  nature  of  the  nonlinear 
element,  are: 

1.  for  a  single-valued,  time-invariant  element: 

if  0  <  k  <   °°  ,  then  -  °°<  q   <   » 
if  k  =  =°  ,    then  0  <,  q   <  « 

2.  for  a  nonlinearity  with  passive  hysteresis: 

0   <   k   <   »    dnd  -    »  <  q    <   0 

3.  for  a  nonlinearity  with  active  hysteresis: 

0   <  k   <   »    and  0    £  q   <  « 

4.  for  a  general  nonlinearity  (time-varying,  and  possibly  with 
hysteresis): 

0   <  k    <   »     and  q  =  0 
If  k  =  »  in  cases  (1),   (2),  or  (4)  above,  then  an  additional  requirement 
must  be  met.    For  all  t  and  for  every  value   em,  there  must  exist  a  finite 
value  un  such  that 

I  u(t)|   su„<Mf     |e(t)|  <  eB  (12) 

It  can  be  noted  immediately  that  the  term  "absolutely  stable"  has 
been  replaced  by  "absolutely  control  and  output  asymptotic".    This  is 

16 


mathematically  the  least  restrictive  practical  stability  definition.    It 
simply  means  that  the  control  signal  or  the  output  signal  eventually  go 
to  zero. 

A  valuable  engineering  tool  provided  by  this  theorem  is  the  trade-off 
available  between  the  requirements  on  the  linear  element  and  those  on 
the  nonlinear  element.     For  example,  if  the  nonlinearity  is  single-valued 
and  time  invariant,  a  non-zero  q  may  be  chosen,  thus  providing  a  greater 
range  for  G(jw)  • 

The  Popov  criterion  provides  another  valuable  asset  to  the  engineer. 
A  measure  of  the  rate  of  system  damping  is  easily  incorporated  in  the 
basic  theorem.    If  the  linear  element  of  the  system  of  Figure  1  is  output 
stable  of  degree  a  ,  this  simply  means  that,    for  a  >  0  ,    the  system 

output,  in  response  to  either  an  initial  condition  or  an  impulse,  tends  to 

-at 
zero  faster  than  the  function   e        .     For  a<0,  it  means  that  the  output 

at 
multiplied  by  e        ultimately  goes  to  zero. 

Now  it  can  be  said  that  if  the  system  of  Figure  1  is  control  asympto- 
tic of  degree  a.  ,  there  exists  a  real  number  a  such  that  for  every  set  of 
initial  conditions 

00 

(eatu(t)  )2  dt   <  »  .  (13) 

Also,  it  is  output  asymptotic  of  degree  a  if  there  exists  a  real  number  a 
such  that  for  every  set  of  initial  conditions 


00 


(eate(t)  )2dt<  oo    .  (14) 
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In  order  for  the  Popov  inequality  to  reflect  a  measure  of  damping,  the 
linear  element  must  be  output  stable  of  degree  a  ,    and  G(jco-a)  must  be 
substituted  for  G(jtc).    Hence,  for  the  system  of  Figure  1  to  be  control 
and  output  asymptotic  of  degree  a  ,  the  following  inequality  must  be 
satisfied: 

Re  (  (l+jqw)G(jW-a)  )  +  1/k  ^  6  >  0  (15) 

Methods  of  system  design  based  on  the  Popov  criterion  have  been 
developed  by  several  authors.    See,  for  example,   [2]  ,   [3]  ,  and  [4]  . 
These  methods  are  all  relatively  straightforward  and  easily  applied  for 
the  case  q  =  0,  beacuse  of  the  invariance  of  the  Popov  line  with  frequency 
on  the  Nyquist  plot.     However,  for  the  case  q  ^  0,  the  design  of  lag-lead 
compensating  networks  involves  considerable  design  effort,  much  of  it 
devoted  to  trial -and -error  computation. 

A  great  deal  of  this  labor  is  eliminated  in  the  new  design  method 
proposed  here  for  the  large  group  of  systems  in  which  design  requirements 
dictate  using  the  case  q  /0.     In  order  to  retain  the  unparalleled  utility 
of  the  polar-plot  representation,  the  lag-lead  compensating  networks 
are  designed  on  the    G*0o.v)  plane. 
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II.     DESIGN  METHOD 

A  sample  design  objective  is  first  illustrated  in  terms  of  a  particular 
system.    Then  the  general  design  method  is  developed. 

Assume  that  the  linear  element  of  Figure  1  is  output  stable  and  has 
the  Nyquist  plot  shown  in  Figure  4.     Suppose  also  that  the  system  design 
requires  the  nonlinear  element  to  be  a  saturating  amplifier  with  the  input- 
output  characteristic  of  Figure  5  . 

The  generalized  Popov  criterion  for  a  general  nonlinearity  (q=0)  gives 
a  maximum  value  of  k  =  1.67  for  the  upper  bound  on  the  Popov  sector  to 
guarantee  control  and  output  asymptoticity  of  the  system.    This  is  shown 
in  Figure  4.    However,  since  the  particular  nonlinearity  of  this  system 
falls  under  Case  1,  a  non-zero    q    may  be  chosen  to  obtain  a  less  con- 
servative result.    In  order  to  do  this,  the  simple  transformation  is  made 
from  the  G(jw)  plane  to  the  G*(jk>)  plane  by  multiplying  the  ordinates  of 
the  G(ja')  curve  by  their  respective  frequency  values.     The  result  is  the 
modified  frequency  response,  G*(joe).     Assume  that  G*(jk')  has  the  form 
shown  in  Figure  6.     Drawing  a  tangent  line  to  the  G*(j&v)    curve  (see 
Figure  6)  at  the  intersection  of  the  G*(j&<)  curve  and  the  negative  real 
axis  gives  a  value  for  q  of  0.1  and  a  value  for  k  of  2.0. 

The  system  nonlinearity  (see  Figure  5),  however,  requires  a  Popov 
sector  with  an  upper  bound  of  at  least  9.0  in  order  to  guarantee  a  stable 
system  with  the  chosen  amplifier.     Hence  the  system  must  be  compen- 
sated.   Excessive  compensation  and  a  too  conservative  solution  can  be 

avoided  by  choosing  a  non-zero  q.    The  design  of  lag  or  lead  cascade 
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ImG*0'u;) 


ReG*(jw) 


-1/k  =  -0.5;  k  =  2.0 


Fig .  6 .    Polar  plot  of  the  linear  element  modified  frequency 
response. 
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compensation  of  the  G*  plane  for  systems  of  this  type  can  be  accomplished 
by  the  method  described  below  with  an  appreciable  decrease  in  labor  over 
previous  methods.    The  compensator  is  to  be  inserted  in  the  output  of  the 
nonlinear  element  in  the  system  of  Figure  1.    Before  proceeding  with  the 
compensator  design,  several  preliminary  relationships  which  are  instru- 
mental to  the  design  method  are  given. 

One  of  the  basic  difficulties  of  designing  in  the  G*(j<*0  domain  is  the 
fact  that 

[G(jo))Gc(ja))]*  ^G*(jto)Gc*(ju>)  (16) 

where  Gc(jw)   is  the  transfer  function  of  the  compensator.    In  order  to 
show  later  how  this  difficulty  may  be  minimized,  the  exact  relationship 
is  derived. 

G(jw)  =  ReG(jw)  +  jlmG(joj) 

Gc(joj)  =  ReGcG6i>)  +  jImGc(jOJ) 

(17) 
G*0o>)  =  ReG(jo))  +  jO)ImGGa.O 

Gc*0w)  =  ReGcGcu)  +  jOJlmGcOw) 

Multiplying  GGw)GcGw)/ 


GGw)GcGw)  =  ReG(jto)ReGc(joj)  -  ImG0w)ImGc0a0 
+  j[ReGGo>)lmGcGo>)  +  ReGcG  w)lmGG  0))] 
/.   [G(jw)Gc(jO))]  *  =  ReG(ja))ReGc(ju>)  -  ImGGu>)ImGcGw) 
+  jw(KeG(jaj)ImGcGw)  +  ReG0G w)ImG(j u) ]  . 
Multiplying  G*G  w)Gc*(jw) , 

G*Gw)G0*(jaO  =  ReGOo>)ReGc(jw)  -  W3ImGGu))ImGcGuj) 
+  ja>[ReG(ju>)ImGc(ja>)  +  ReG0(jw)ImGGa))]  . 
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(18) 


(19) 


(20) 


Comparing  (19)  and  (2  0), 

[G(jw)Gc(ju>)]  *  =  G*(jw)Gc*(J05)  +  (aJ2-l)ImG(jw)ImG0(j^)   .        (21) 
The  frequency-response  plots  of  the  compensation  networks  are 
essential  to  the  design  procedure.    The  transfer  function  of  a  lag  compen- 
sation network  is  given  by 

Gc(s)  =    YT^2     where    Tl>Ts  (22) 

and  that  of  a  lead  network  is  given  by 

Gc(s)  =    ]  +STl      where    Tl   >    Ts  .  (23) 

1  +  st2 

The  Nyquist  plots  of  the  two  networks  are  shown  in  Figure  7. 

The  geometric  configuration  of  the  modified  frequency  responses, 

Gc*(joO,  is  now  derived  for  the  lag  filter  only.    The  lead  filter  gives  the 

same  result  only  with  ra   and   r2  reversed.     For  the  lag  filter: 

Gc*(jW)  =      1+t^IP    +     J         "3<T%"  'i      ■  (24) 

1   +  0?    Ti  1   +  Ui    Ti2 

•    X  =     )  +  w;T>;»     ;       Y    -      f  (T%-  % »     .  (25) 

1    +  tt'2  T*  1    +  03    Tx2 

Solving  for  u>2  in  the  X-expression: 
(1  +  w3  Tx2  )X  =    1  +  t»2  Ti  T2 
X  +  a)3  tx2  X  =    1  +  W2TiT2 
0)2Ti2X  -4JaTi  TS  =    1    -  X 

.-.  **  =  t4itxt  (26) 
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Substituting  into  the  Y-expression: 
1  -  X 


Ti    X  -   TXTZ 
Y=    " 


*     (T3    -    1"!  ) 


.  1    ~  X 2 

1  +   -ax  —    , V 

'  1     A         Tl    T2 


(i-x)  (TP-  rj 
Ti2X  -  tx  r2 


T£X  -  txt3  +  (l-X)iy 

Ti2X  -   Tl   T2 


(i  -x)  (t2  -  n ) 


n2X  -  Tl  T3  -  ra2X  +  Tl2 


(1-X)   (Ta-  Tj 
Ti    (TX   -   T2  ) 

(X-l)   (Ta-Tx) 
Ti    (Ts-Ti) 

Therefore , 

Y=  1/tj  (X-l)   .  (27) 

Similarly  for  the  lead  filter, 

Y=  l/r2  (X-l)   ,  (28) 

which  is,  of  course,  the  equation  of  a  straight  line.    This  very  conven- 
ient result  is  instrumental  to  the  design  method.    In  order  to  compute  the 
end  points  of  the  line  (24)  is  simply  evaluated  at  0)=  0  and  oo=  °°  .    The 
results  for  both  filters  are  shown  in  Figure  8.    The  frequencies    a?  =  1/tx 
and  co  =  1/t2  occur  at  the  geometric  means  of  the  respective  polar  plots. 

The  ideas  developed  so  far  can  now  be  applied  to  the  system  with 
the  modified  frequency  response  shown  in  Figure  6.     Suppose  we 
decide  to  compensate  this  system  with  a  lag  compensator  to  force  the 
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compensated  locus,     [GO  w)Gc(jO})]  *,  to  intersect  the  negative  real 
axis  at  a  point    ^  -  0.1  and  at  a  certain  crossover  frequency,     wc   . 
This  would  give  a  Popov  sector  of  at  least    (0,   10)  which  would  guar- 
antee a  stable  system  with  the  required  amplifier,  provided  that    a 
Popov  line  of  non-zero  q  can  be  drawn  through  the  point    (-0.  1,0). 

The  procedure  to  accomplish  this  alteration  of  the    G*(jkv)  locus 
is  first  outlined  briefly,  then  in  more  detail.    The  desired  negative 
real  axis  crossover  point  for  the  compensated  locus  is  designated  as 
P0  .    Then  the    l/G*Gw)    and  the    G0*(joo)/P0    loci  are  constructed  so 
as  to  intersect  at  the  point    Pa    where  the  value  of  the     1/G*(jw) 
locus  corresponds  to  the  chosen  crossover  frequency,    w0/  and  the 
parameters,     tx  and   t2    ,    are  then  computed  so  that  the  frequencies  of 
the  two  curves  are  equal  at  the  intersection  point.    This  construction 
designs  the  compensator  since  it  guarantees  that 

1/G*(jw)  =Gc*Gw)/P0      at    w  =    Wo 
or 

G*(ja>)G0*(jw)  =  P0      at    u>  =   wc 

where  G*(jw)G0*(jw)    and  [Gc(j w)G(j w)l  *   are  related  by  (21).    This 
relationship  is  discussed  in  detail  later. 

To  show  how  easily  this  is  accomplished  in  practice,  the  example 
of  this  section  is  used.     Let  P0  =   -0.1  and  coc  =   1.0.    The  first  step 
would  be  to  construct  the  l/G*G'w)  curve.     However,  this  is  unnecessary 
since  only  one  point  of  the  1/G*(jw)  curve  is  actually  used,  the  one 
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corresponding  to    03c  =   1.0.    The  coordinates  of  this  point,   Px,  are  easily 
computed  and  plotted  by  inverting  the  complex  number  read  off  the  G*(jU>) 
locus  for  ccc  =   1.0. 

The  next  step  is  to  construct  the  Gc*(jto)/P0  plot  by  simply  rotating 
the  Gc*^^)  plot  by  180°  and  placing  the  w  =   0  point  of  the  compensator 
at  1/P„ .    The  slope  of  this  straight-line  plot  is  adjusted  so  that  the  line 
passes  through  the  point  P± .    Thus  Px   corresponds  to  the  intersection  of 
the  1/G*(oj)  and  the  G£(ju>)/P0    loci,  and  the  frequency  at  this  point 
must  be  the  same,  i.e.,     o)=  (jici  on  both  curves.     Refer  to  Figure  9. 

Equation  (27)  gives  the  equation  of  the  lag-compensator  plot.    By 
comparing  this  with  the  standard  point-slope  form  of  a  straight  line,  it 
can  be  seen  that  the  slope  is  1/tx  .    Hence    rx    can  be  determined  by 

■ 

simply  computing  the  slope  of  the  line  graphically.     Then    t2  is 
calculated  from  either  part  of  (25)  using    u$  =   wc.    The  coordinates,    X 
and  Y,    are  simply  given  by  the  real  and  imaginary  parts  of    Px  /P^    . 
This  is  a  very  easy  calculation  since  the  value  of  the  squared  parameter 
(in  this  case,     tx  )  is  always  substituted  numerically,  whether  the 
filter  is  lag  or  lead.     Using    w  =  W0  in  the  calculation  of    t2  guarantees 
that  the  frequency  of  the  G0*(ju>)/P0  function  at  Px  is  the  same  as  that 
of  1/G*(jw)   at  Pj_  . 

To  avoid  a  sign  error  when  computing  Px/P0,  it  should  be  noted  that 
X  and  Y  are  always  positive  for  a  lead  compensator  and  X  is  positive  and 
Y  is  negative  for  a  lag  compensator.    It  should  also  be  noted  that  there 
are  slope  limitations  on  the  straight-line  compensator  plots.    Infinite 
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slope  corresponds  to  the  pole  and  zero  which  cancel  each  other  while 
zero  slope  corresponds  to  infinite  pole-zero  separation. 

The  frequency  response  curve  of  interest  in  the  system  design  is  the 
[G(ja)Gc(j6o)  ]  *  curve.    The  plot  given  by  the  compensation  method 
explained  above  gives  the  G*{]Ui)Gc*(jUi)  curve.    Equation  (21)  relates 
these  two  quantities.    Note  that  since  toc  was  chosen  to  equal  1.0  in  the 
above  description,  the  two  curves  cross  the  negative  real  axis  at  the 
same  point.    This  is  not  the  case  in  general. 

Assume  that  the  two  resulting  curves  are  as  shown  in  Figure  10. 
Here  a  non-zero  q  line  can  be  drawn  through  the  point  (-0.1,0)  and  the 
system  design  is  completed.    In  an  actual  design  problem  there  may  be 
other  performance  specifications,  such  as  minimum  bandwidth,  which 
must  be  met  before  the  design  is  complete.    These  additional  requirements 
will  influence  the  choice  of  toc  and  P0 

In  order  to  generalize  the  method  to  any  0Je/   (21)  is  repeated  for 
reference. 

[G(jw)Gc(jw)]*  =  G*(j^)Gc*(jw)  +  (a>2-l)ImG(jO>)ImGc(jC0).  (21) 

For  simplicity,  let  (<A)2-l)lmG(jto)ImGc(juO    -  T.    It  can  be  noted 
immediately  that 

T 

)  when  in         — |  0  and  »    since  ImG0(jw)   =   0  (2  9) 

value  for  which  ImG(j&0   =   0  . 
For  the  compensated  system  to  be  worsened  by  the  effect  of  T,   (i.e. 
the  [G(ju>)Gc(jaj)3  *  curve  be  more  negative  than  the  G*(j to)Gc*(jw)  curve 
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Popov 
line 


[G(jw)GcGw)3* 


\ 


Im[G0oi)Gc0uj)]* 


G*(jO>)Gc*(jw) 


0J  = 


Re[G(jw(Gc0w)]* 


Fig.  10.    Modified  frequency-response  locus  for 
compensated  system. 
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at  the  negative  real  axis  crossover),  the  term  T  must  be  negative  for 
W=   Wc  . 

For  a  lag  network,  since  ImGc(jOo)  equals  zero  or  a  negative  value  for 
all  w,  this  can  occur  only  under  the  following  two  conditions: 
if    60c    <  1  and  ImG(j60o)   <  0 
or  if    60c    >  1  and  ImG(jwc)    >0.  (30) 

For  a  lead  network,  since  ImG0(j6o)  equals  zero  or  a  positive  value  for 
all  co,    T  can  be  negative  only  under  the  following  conditions: 
if  60c  <    1  and  ImG(j60c)  >  0 
orif  wc>    1  and  ImGOtOc)    <0  .  (31) 

Using  the  methods  outlined  above,  a  rapid  design  procedure  can  now 
be  stated.    Since  the  subject  of  compensation  is  treated  in  most  control 
system  design  textbooks,  no  attempt  will  be  made  in  this  paper  to  sum- 
marize the  various  helpful  hints  given  by  other  authors  for  compensation 
design  of  linear  systems  using  lag  and  lead  networks.     Since  the  design 
method  outlined  here  involves  reshaping  a  linear-system  frequency- 
response  curve,  it  is  assumed  that  the  designer  has  a  feel  for  the  steady- 
state  and  transient  effects  of  compensation  on  linear  system  performance 
and  is  cognizant  of  the  advantages  and  disadvantages  of  lag  and  lead 
compensators.    See,  for  example,   (5).    With  this  in  mind,  the  design 
steps  may  be  listed  as  follows: 

1.     Plot  the  polar  curves  for  G(jco)  and  G*(jto)  as  60  varies  from  0 
to  °° .    The  G*(ju>)  plot  is  easily  made  by  multiplying  the 
ordinates  of  the  G(j6o)  locus  by  their  respective  frequencies. 
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2.  From  the  required  design  specifications  and  system  compo- 
nents choose  the  type  of  compensation  (lag  or  lead)  and  Wc  . 

3.  Use  (30)  or  (31)  and  the  design  specifications  for  the  initial 
choice  of  P0  . 

4.  Design  the  compensator. 

5.  Sketch  the  G*(jO>)Gc*(jw)  locus  by  evaluating  the  compensator 
of  Step  4  for  a  minimum  number  of  frequencies.     Sketch,  if 
necessary,  the  [G(ja.0Gc  (ja*)]  *  locus  using  (21). 

6.  Evaluate  the  design  in  view  of  the  specified  performance 
criteria. 

If  multiple  sections  are  required  the  process  is  simply  repeated 
using  the  [G(jaOGc(jo.')  ]  *  curve  in  Step  1.     Since  the  design  method  is 
virtually  identical  for  lead  and  lag  compensators,  a  detailed  discussion 
of  lead  network  design  is  left  to  the  examples  section.     It  should  also  be 
noted  that  the  design  method  introduced  above  is  quite  flexible.     No 
additional  complexity  is  introduced  when  a  P0  is  desired  off  the  negative 
real  axis.     l/Po  is  simply  designated  at  the  desired  point  and  the 
Gc*(ju>)  line  is  rotated  the  required  number  of  degrees. 


33 


III.     EXAMPLES 

To  illustrate  how  the  preceding  results  might  be  used,  the  following 
examples  are  given. 

Example  1 .    It  is  desired  to  compensate  the  system  shown  in  Figure 
11  so  that  the  output  signal  will  be  asymptotically  stable  of  degree  0.1. 

To  simply  guarantee  output  asymptoticity  with  the  specified  type  of 
nonlinearity,  the  minimum  Popov  sector  must  be  [0,5]  .    However,  for 
the  system  to  be  asymptotically  stable  of  degree  0.1,  the  Popov  line 
(for  some  q)  must  intersect  the  negative  real  axis  at  a  value    ^-0.2  and 
must  not  be  closer  than  0.1  horizontal  units  to  the  G*(jco)  curve  at  the 
closest  point. 

The  G(joo)  and  G*(jui)  curves  are  shown  in  Figures  12  and  13.    A  lag- 
type  compensator  and  an  wc  =1.0  are  selected.    To  avoid  the  possibility 
of  a  borderline  design,  a  P0  of  -0.08  is  used,  as  opposed  to  0.1.    The 
graphical  design  of  the  compensator  is  shown  in  Figure  14.    The  1/P0 
point  is  at    1/-0  .08  =  -12  .5  .    To  locate  point  Px  =  1/G*(jl),  note  on 
Figure  13  that  G*(jl)  =  -0.4-J0.2  from  which  1/G*(jl)  =  -2+jl  as  shown. 
The  parameter  values,  \  and  t2  ,  are  computed  using  the  slope  of  the 
Gc*(juO/P0   Plot  (see  Figure  14)  and  Equation  (25).    The  resulting  compen- 
sator is  given  by 

G<(s)  ■  frr^i  •  <32> 

Using  a  minimum  number  of  points  the  G*(j0t>)Gc*(jo>)  locus   (see 
Figure  15)  is  easily  sketched.    Then,  using  (21),  with  (30)  as  a  check, 
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ImGOoo) 


EeG  (j  u>) 


-0.1 


-0.2 


-  -0.3 


Fig.   12.    G(jtii)  locus  for  Example  1 
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Fig.   13.    G*(joj)  locus  for  Example  1 
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u>  =  2 


G*Gu>)Gc*(joj)   plane 


Fig.  15.    G*(j60)Gc*Go>)   locus  for  Example  1 
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the  [G(jw)Gc(jto)]  *    locus  can  be  quickly  determined.    In  (21),  the 
(a)2-l)ImG0w)ImGc(j6u)  term  can  be  evaluated  by  simply  picking  the  values 
of  ImGOw)  off  the  Nyquist  plot  of  G(jw)  and  using  the  expression  for  Y 
in  (2  5)  divided  by    to  for  Im  Gc(j<o). 

In  this  problem,  the  G*(jto)Gc*(ju))  and  the  [G(juOGc(jto)]  *  loci  have 
the  same  negative  real  axis  crossover  at    t*>c  =  1  and  are  virtually  identi- 
cal at  other  points.    The  Popov  line  is  as  shown  in  Figure  15  guaranteeing 
control  and  output  asymptoticity  of  degree  0.12. 

Example  2  .    Consider  the  nonlinear  control  system  with  time  delay 
shown  in  Figure  16.    Hsu  and  Meyer  [2]  studied  this  system  and  obtained 
the  following  results: 

1.  For  a  class  2  or  class  4  nonlinearity  (q  =  0) ,  the  Popov  criter- 
ion is  satisfied  for  0  <k  <  1.85  (see  Figure  17). 

2.  For  a  class  1  or  class  3  nonlinearity  (q  ^  0) ,  the  Popov  criter- 
ion is  satisfied  for  0  <  k  <  1.98  (see  Figure  17). 

Now  it  is  desired  to  compensate  this  system  to  permit  a  Popov  sector  of 

[0,5]   for  a  class  1  or  class  3  nonlinearity. 

The  G*(jto)  locus  is  shown  in  Figure  17.     Lag  compensation  is  chosen 

with  an   u>c  =    1.8.    Equation  (30)  quickly  shows  that  if  a  P0  is  chosen 

which  gives  a  minimum  Popov  sector,  the  initial  design  will  not  produce 

an  unstable  system  but  will  actually  give  a  greater  degree  of  stability. 

Hence  a  P0  of  -0.2  is  used.     Designing  the  compensator  as  before  (see 

Figure  18)  gives  the  result 
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Fig.   17.    G*(j(jo)  locus  for  Example  2. 
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Therefore: 

1.  For  a  class  2  or  class  4  nonlinearity  (q  =  0) ,  the  Popov  criter- 
ion is  satisfied  for  0  <k  ^  4.45   (see  Figure  19).    This  example 
illustrates,  as  a  sidelight,  how  easily  the  design  method  may 
also  be  applied  to  the  case  when    q  =  0,  although  the  advan- 
tages are  not  quite  as  great  as  when  q  ^0. 

2.  For  a  class  1  or  class  3  nonlinearity  (q  /  0) ,  the  Popov  criter- 
ion is  satisfied  for  0  <  k    ^5  (see  Figure  19.) 

The  permissible  Popov  sector  will  actually  be  slightly  greater  than  above 

as  indicated  by  (30).    If  exact  sector  limitations  are  desired,  these  are 

quickly  obtained  by  using  (21)  and  (30)  to  obtain  the  [G(j td)Gc(jo>)]  * 

locus . 

Example  3  .    Consider  the  actuator  or  indirect  control  system  of 
Figure  20.     Hsu  and  Meyer  [2]  also  studied  this  sytem  using  the  gener- 
alized Popov  theorem.    Their  intent  was  to  determine  the  critical  feed- 
back gain  F  =  Fc  such  that  for  feedback  gains  F  ^  F0,  the  output  signal 
y(t)  would  asymptotically  approach  a  steady-state  value  dependent  on 
the  constant  input  value,  r. 

The  modified  block  diagram  of  this  system  is  shown  in  Figure  21  . 
This  zero-input   model  is  possible  because  the  constant  input,     r,     can 
be  combined  with  the  constant-initial-condition-response  term  with  no 
loss  of  generality . 

Using  the  modified  frequency  response  function, 

G*(j<j)  =  G2*(j<*>)  -  JF  (34) 

where 

r*(<,A-  24  +  g)2  36  -  gj5    „___ 

G2"VJU/"    "  36  +  i3w2+  w4  J      36  +  i3u.'2  +  o>4      '  ^' 
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For  asymptotic  stability,  the  Popov  criterion  with   k.  =  »  then  requires 
that 

ReG*(ja«)-qImG*(jcl>)  =  ReG2*G^)-q(ImG2*(Jkl)-F)  >  0. 

G2*(JtL')  is  plotted  in  Figure  22,     Since  the  ImG2*(Jw)  has  a  maximum 
value  of  approximately  0.006,      the  Popov  criterion  is  satisfied  with  some 
q  >0  if  F  >  0.006  and  the  system  is  then  control  and  output  asymptotic. 
Hence  e(°°)    =   0   and 

y(oo)  =   f+Gi(o)    "   TT7~  r  (36) 

for  this  system.     Therefore  the  range  of    F    for  which  the  system  is 
control  and  output  asymptotic  is  extremely  important  since  F  directly 
affects  the  steady-state  error. 

Digital  simulation  of  the  system,  with  F  =  0.006,   provided  results 
agreeing  with  the  theoretical  values.     The  response  of  the  simulated 
system  to  an  input  value  of  r  =  i  .  0  is  shown  in  Figure  23  . 

Assume  now  that  design  requirements  specify  a  steady-state  error 
which  requires  the  system  of  Figure  20  to  be  asymptotically  stable  for 
F  =  0.004.     Hence  the  system  must  be  compensated  so  that  the  maximum 
value  of  the  ImG2*(ja.O    is  less  than  0.004.     For  comparison  purposes, 
the  unit  step  response  of  the  uncompensated  system  with  F  =  0.004  is 
shown  in  Figure  24.     Lead  compensation,    o)c  =   8.8,  and  P0  =  -  .017  are 
chosen.     The  compensator  design  is  shown  in  Figure  25  with  the  result 

n   (a\  l  +  °-ls  /„\ 

Gc(5)    "     1+  0.079s  (37) 
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Fig.  23.    Response  of  the  system  of  Figure  20  to  a  unit  step 
input.    F  =  0.006. 
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Fig.  24.    Response  of  the  system  of  Figure  20  to  a  unit  step 
input.    F  =  0.004. 


51 


60  -40 

1/Pn  =   "59 


Gc*(jw)/Pc 


■20 
_i — 


ImGc*(ju«)/Pc 


-ReGc*(ja')/Pc 


Fig.  25.     Lead  compensation  design  plot 
for  Example  3 . 
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The    Gc*(ja>)G2*(jw)    locus  is  shown  in  Figure  26. 

It  should  be  noted  that,  in  this  case,   the  T  term  of  (29)  is  of  no 
consequence  since  it  is  only  real  and  will  have  no  effect  on  the  maximum 
imaginary  value  of  the  compensated  curve.    Hence 

maxim  [Gc(j  w)G2(jo))3  *   =  maxim   G0*(Jtt)Gs*(j  U))  (38) 

and  the  design  is  completed.     The  response  of  the  compensated  system 
with  F  =  0.004  and  r  =  i  .  0  is  shown  in  Figure  2  7  . 
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Fig.  26.    G2*(jw)Gc*(jco)  locus  for  Example  3. 
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Fig.  27.    Response  of  the  compensated  system  of  Example  3 
to  a  unit  step  input.    F   =   0.004. 
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IV.     CONCLUSION : 

A  design  method,  based  on  the  Popov  stability  criterion,  has  been 
developed  for  the  large  class  of  nonlinear  systems  which  can  be  modeled 
in  the  configuration  of  Figure  1.    The  design  procedure  uses  the  modified 
frequency-response  polar  loci  for  the  linear  portion  of  the  system.    This 
permits  the  engineer  to  utilize  the  clear  physical  picture  inherent  in  the 
Nyquist  plot. 

As  a  major  step  in  the  development  of  the  design  procedure,  it  was 
shown  that  the  modified  frequency-response  polar  plots  for  lag  and  lead 
compensation  networks  are  straight  lines  .    This  provided  the  basis  for 
the  development  of  a  rapid,  straightforward  design  method  which  could 
be  quite  useful  in  engineering  practice. 
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